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ABSTRACT 

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer 
simulations for supernovae (SNe) of type lib, lb, and Ic that result from the terminal explosion 
of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic en¬ 
ergy of w 1.2x 10 51 erg, but characterised by different ejecta masses (2-5 M 0 ), composition, 
and chemical mixing. The type Ilb/Ib models correspond to the progenitors that have retained 
their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has 
lost its helium shell, but retains 0.32 Mq of helium in a CO-rich core of 5.11 Mq . We discuss 
their photometric and spectroscopic properties during the first 2-3 months after explosion, and 
connect these to their progenitor and ejecta properties including chemical stratification. For 
these three models, Arnett’s rule overestimates the 56 Ni mass by « 50% while the procedure 
of Katz et al., based on an energy argument, yields a more reliable estimate. The presence 
of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was 
under-abundant in helium. As noted by others, the 1.08/rm feature is a complex blend of C I, 
Mg II, and He I lines, which makes the identification of He uncertain in SNe Ibc unless other 
He I lines can be identified. Our models show little scatter in (V — R) colour 10 d after /Aband 
maximum. We also address a number of radiative transfer properties of SNe Ibc, including the 
notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum for¬ 
mation, blackbody fits and “correction factors”. 

Key words: radiation hydrodynamics - radiative transfer - stars: supernovae - stars: evolu¬ 
tion - stars: binaries 


1 INTRODUCTION 

SNe Ilb/Ib/Ic are understood to arise from the collapse of the iron 
core of massive stars having lost most (type lib) or all (types 
lb and Ic) of their envelope exterior to the He core. Wolf-Rayet 
stars are therefore the natural progenitors of these transient phe¬ 
nomena. However in the 80’s, after the first well sampled light 
curves of type I core-collapse SNe we re obtained, it became evi- 
dent thaMhejr e jecta were of low mass fensm an & Woqsjev 11 98i~j ; 
Wooslev et al.ll 19951) . in accord with the binary evolution scenario 
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and moderate progenitor masses IWelistem & Lanuen 11 99SI) . As 
wind mas s loss rates for Wolf-Rayet stars were revised downwards 
(see, e.g., lHillierlfT99ll . iNugis & Laniers 2000h and inferences for 


'essart et 


ah] 


201l| J Drout et al. 2011; Bersten et al. 2012 

: Hachinger et ah 2012]; 

Benvenuto et alj|20l5: Ergon et al.1 20141: 

Jerkstrand et ah 20151), 


the binary star evolution channel became the ma in scenario for the 
production of type I core collapse SNe (see, e.g., lYooijl2015t) . 


An extensive study of binary evolution, and its im plications 
for SN types, was made by IPodsiadlowski et ahl dl992lh As part 
of their study, they estimated probabilities for various evolution¬ 
ary scenarios for massive stars with masses between 8 and 20 Mq. 
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2 Luc Dessart et al. 


Table 1. Summary of progenitor properties at the time of core collapse. The 
last column gives the iron-core mass. 


Model 

M; 

[M 0 ] 

M f 

[M 0 ] 

L* 

[erg s^ 1 ] 

T e ff 

[K] 

R* 

[cm] 

Mco 

[M 0 ] 

MFe 

[Mq] 

3p65 

16.0 

3.65 

2.16(38) 

21700 

1.17(12) 

2.03 

1.43 

5pll 

60.0 

5.11 

5.42(38) 

130 000 

5.19(10) 

3.64 

1.52 

6p5 

25.0 

6.50 

6.40(38) 

68 700 

2.01(11) 

4.21 

1.48 


They concluded that a complete understanding of massive star evo¬ 
lution will only be obtained when we understand binary star evo¬ 
lution. This also requires a ce nsus of binaries and their parameters. 
A more recent study is that of lClaevs et alj J201 ll) who studied the 
roles of binaries in producing type lib SNe. 

In this first paper of a series we perform 1-D non-Local- 
Thermodynamic-Equilibrium (non-LTE) time-dependent radiative- 
transfer simulations for supernovae of type lib, lb, and Ic that re¬ 
sult from the terminal explosion of the primary star in a close- 
binary system. The progenitor models are taken from the study 
of lYoon et al.l J2010I ). which treats the evolution from the main 
sequence until neon burning. The later stages of evolution, and 
the subseque nt explos ion is don e with KEPLER (Woosley et al., in 
preparation; W eaver et all 19781 : Section©. 

_ This e n deavour e x tends beyond the previous studies of 

iDessart et alj J201 1 LI 2 OT 2 I) by using better physics, a broader range 
of progenit ors, and a broader range of explosion energies and 56 Ni 
mass. In IDessart et alj d201 lh. non-the r mal e ffects and chemical 
mixing were ignored. In Pessart et alj (1 201 2l) . the focus was on 
mixing alone, based on only two explosion models with the same 
energy. Both progenitor models had a large helium mass. Here, 
the full set comprises about 30 model sequences covering a much 
broader parameter space. We compute full time sequences in non- 
LTE, with allowance for time dependence, non-thermal effects and 
non-local energy deposition (Section©. Throughout this series, we 
limit the discussion to the photospheric phase. The nebular phase 
properties, because they are so strongly connected to the nature of 
mixing, will be treated in a separate study. Our goal is to charac¬ 
terise the light curve and spectral properties of SNe Ilb/Ib/Ic and 
correlate them to the explosion model and to the progenitor star. 
This includes a discussion of the bolometric light curves, multi¬ 
band light curves, colour evolution (Section 0, spectral evolution 
(Section©, the properties of near-IR He I lines (Section©, the in¬ 
fluence of mixing (Section©, the properties of spectrum formation 
and the notion of a photosphere in type I SNe (Section©, the in¬ 
ference of the ejecta expansion rate and the constraint on chemical 
stratification (Section©, and the correction factors to be used when 
approximating the synthetic spectra of SNe Ilb/Ib/Ic near maxi¬ 
mum light with a blackbody (Section©. 

In Section [TO] we summarise the results of the present paper. 
In Paper II, we will do a statistical study based on the entire grid 
of about 30 model sequences to examine correlations/distinction 
between the models. In Paper III, we will use our grid of models to 
compare to observations of SNe Ilb/Ib/Ic. 


Table 2. Summary of ejecta properties. The last three columns give the 
ejecta velocity that bounds 99% of the corresponding species total mass. 
The integration is done inwards in velocity space for H and He, and out¬ 
wards for 56 Ni. 


Model 

M r 

[M 0 ] 

Me 

[M 0 ] 

f^kin 

[B] 

Vg9,H 

[km s _1 ] 

Vg9,He 
[km s -1 ] 

V99,Ni 

[kms" 1 ] 

3p65Axl 

1.43 

2.22 

1.24 

1.34(4) 

3.06(3) 

7.00(3) 

3p65Ax2 

1.43 

2.22 

1.22 

1.24(4) 

2.62(3) 

1.01(4) 

5pl lAxl 

1.57 

3.54 

1.25 


1.68(3) 

6.20(3) 

5pl 1 Ax2 

1.52 

3.59 

1.29 


2.71(3) 

9.18(3) 

6p5Axl 

1.53 

4.97 

1.26 


3.28(3) 

5.98(3) 

6p5Ax2 

1.55 

4.95 

1.25 


2.55(3) 

8.60(3) 


until the end of neon burning, including differential rotation, tides 
and mass and angular momentum transfer, for a limited set of 
close-binary systems (orbital period of «4d), covering a range of 
primary-star masses (12 to 60 Mq) and mass ratios (>1-1.5). These 
simulations were then evolved to core collapse (Woosley et al., in 
preparation), and subsequently exploded with KEPLER using a pis¬ 
ton. 

We have performed a grid of explosion models, cor¬ 
responding to the binary models in the Yoon et al. study 
that are numbered Yoon_13 (the final mass of the primary 
star, Mi,f, is 3.0 Mq), Yoon_9 (Mi,f = 3.65M 0 ), Yoon_35 
{Mi s =4.64M 0 ; Z =0.004), Yoon_3lQ (M M =5.11 M 0 ), and 
Yoon_30 (Mi,f =6.5M 0 ). All models use the solar composition 
initially, except model Yoon_35. This choice includes models with 
residual surface hydrogen (Yoon_13, 9, 35), no hydrogen but high 
helium and nitrogen abundances (Yoon_30), and little helium but 
high abundance of carbon and oxygen (Yoon_31). 

In this paper, we limit the discussion to just three progenitor 
models, 3p65, 5pll, and 6p5 (Table©. Their diversity in envelope 
composition leads them to produce, after explosion, the three com¬ 
mon type I core-collapse SN lib, lb, and Ic. Models 3p0 and 4p64 
will be discussed in Papers II and III. 

Upon reaching core collapse, each model was exploded with 
KEPLER to produce ejecta with an asymptotic kinetic energy of 0.6 
(suffix C), 1.2 (suffix A), 2.4 (suffix B), and 5x 10 51 erg (suffix D). 
In practice, the lower mass models were limited to lower energy ex¬ 
plosions (0.6-2.4X 10 51 erg). Explosive nucleosynthesis is treated 
in KEPLER and therefore these simulations naturally produce some 
56 Ni. Higher energy explosions exhibit higher post-shock tempera¬ 
tures, which favour 56 Ni production. The full set of explosion mod¬ 
els is characterised by an original 56 Ni mass of 0.05-0.2 M 0 . Here, 
we only discuss “A” models, characterised by a kinetic energy at in¬ 
finity of « 1.2x 10 51 erg. 

These explosion mod els were then remapped into Id dLivnel 
ll993l ; lDessart et al.l2010al lbh at 10 4 s and evolved to 10 5 s. This ad¬ 
ditional step is to explore a different mixing approach from that 
employed in KEPLER. In our 1-D approach, we take the multi- 


dimensional process of mixing dFrvxell et alj 1991: 

Kifonidis et al. 

200d: Hammer et al. 2010i; Wonawathanarat et al.l 

2014) into ac- 


count by applying a simple algorithm. Starting from the ejecta 


2 NUMERICAL SETUP AND INITIAL CONDITIONS 

2.1 Pre-supernova evolution, explosion, and supernova 
radiative transfer 

Yoon et alj d20ld) computed the evolution from the main sequence 


1 This m odel has a final m ass of 5.11 Mq, and not 4.96 Mq as given in 
Table 1 of I Yoon et alj fcOlQh . This is because the remapping of this model 
into KEPLER was done at an earlier stage and the final evolution to core 
collapse neglected wind mass loss. 
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Figure 1. Ejecta composition at « 2.5 d after explosion for models 3p65Axl (top row), 5pl 1 Axl (middle row), and 6p5Axl (bottom row). In the left column, 
we show the mass fraction versus velocity for a selection of species. The dashed line corresponds to 56 Ni. The top axis gives the lagrangian mass coordinate. 
For clarity, the curve for 56 Co is not shown; at 2.5 d it amounts to « 1/4 of the 56 Ni mass. In the right column, we show for each mass shell m the cumulative 
mass fraction Xj(ra), where i covers species H, He, C, N, O, Ne, Mg, Si, S, and Ni (for Ni, we show the current total mass fraction, which therefore 
accounts for decays). We also show the velocity (dashed line, right axis). This illustration reveals more clearly the relative fraction occupied by abundant 
species. The variable size of the mass buffer between the He-rich layers and the 56 Ni rich layers is key to understand the production of He I lines. 
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Table 3. Cumulative yields for our grid of models at ^2d after explosion. For 56 Ni, we give the original mass, i.e., prior to decay. Models with different 
mixing show slightly different levels of fallback, and consequently show slightly different yields. 


Model 

H 

[M 0 ] 

He 

[M 0 ] 

C 

[M 0 ] 

N 

[M 0 ] 

O 

[M 0 ] 

Si 

[M 0 ] 

S 

[M 0 ] 

Ca 

[M 0 ] 

( 56 Ni) 0 

[M 0 ] 

3p65Axl 

4.99(-3) 

1.49(0) 

9.40(-2) 

1.09(-2) 

3.02(-l) 

7.38(-2) 

2.36(-2) 

3.94(-3) 

7.42(-2) 

3p65Ax2 

4.72(-3) 

1.48(0) 

9.37(-2) 

1.08 (-2) 

3.0K-1) 

7.30(-2) 

2.331-2) 

3.91(-3) 

7.66(-2) 

5pl lAxl 

0 

3.15(-1) 

8.921-1) 

0 

1.42(0) 

1.28(-1) 

4.001-2) 

5.98(-3) 

8.94(-2) 

5pl 1 Ax2 

0 

3.26(-l) 

9.231-1) 

0 

1.44(0) 

1.281-1) 

4.041-2) 

6.13(-3) 

9.46(-2) 

6p5Axl 

0 

1.67(0) 

4.13(-1) 

7.59(-3) 

1.57(0) 

2.12(-1) 

9.31 (-2) 

9.30(-3) 

9.90(-2) 

6p5Ax2 

0 

1.66(0) 

4.11 (-1) 

7.52(-3) 

1.57(0) 

2.14(-1) 

9.40(-2) 

9.51(-3) 

1.02(-1) 


Table 4. Mass fractions for some important species in the outermost mass shell for our grid of models. This shell corresponds to the progenitor surface. The 
same outermost-shell composition is obtained for the x2-model counterparts (hence not shown). 


Model 

A h , s 

^He,s 

^C,B 

Tn,s 

A 0 ,b 

A S i,s 

^S,B 

-^Ca,s 

-^Fe,s 

3p65Axl 

1.480(-1) 

8.332(-l) 

1.980(-4) 

1.320(-2) 

4.831(-4) 

7.352(-4) 

3.651(-4) 

6.442(-5) 

9.086(-4) 

5pl lAxl 

0 

3.611(-1) 

5.05H-1) 

0 

1.100(-1) 

7.362(-4) 

3.651(-4) 

6.442(-5) 

9.097(-4) 

6p5Axi 

0 

9.813(-1) 

4.147(-4) 

1.309(-2) 

3.098(-4) 

7.345(-4) 

3.647(-4) 

6.435(-5) 

9.086(-4) 


base, we make homogeneous all mass shells within 1000 (xl) or 
2000 km s~ x (x2), and sweep this mixing-box through all mass 
shells all the way to the outermost ejecta layers@ This mixing 
box is thus applied hundreds of times (i.e., the number of times 
is equal to the number of grid points in the model). Beca use of 
the o v erwhelming evidenc e that SN Ib/c ejecta are mixed JLucvl 
1 199 ll: De ssart et alj|2012h . we do not include unmixed models. 
But because the magnitude of mixing is very uncertain, we use 
two values corresponding to moderate and strong mixing. Our 
approach enforces a macroscopic mixing (large scale mixing of 
species over a wide range of velocities) as well as a microscopic 
mixing (each mass shell is homogeneous with no segregation of 
species at a given velocity). Multi-dimentional simulations of core¬ 
collapse SNe predict stro ng macroscopic mixing, but v ery limited 
microscopic mixing (e.g., Wongwathan arat et Hu I2oi 4 and refer¬ 
ences there in). O bservational ly, this distinction becomes evident at 
nebular times Qjerkstrand et alj|2015i) . Models with different mix¬ 
ing produce ejecta that suffer a slightly different fallback. This 
leads to small variations in yields (Table [3J. 

Our nomenclature is to name each model according to its final 
mass. For example, model 3p65Axl corresponds to model Yoon_9 
which died as a 3.65 M@ star. It was exploded to yield an asymp¬ 
totic kinetic energy of 1.2 x 10 51 erg, and was mixed with a velocity 
width of 1000 km s _1 . 


At 10 s s, all models are remapped onto the grid of the non- 
LTE r adiative transfer code C MFGEN and evolved until nebular 
times dHillier & Dessartll20 1 2t) . The entire ejecta is treated, and 
each mode l is evolved with the same physics . We treat non-thermal 
processes i Dessart et all201 AlLi et all2012l ) and time dependence 
(iDessart & Hilliej2008l) . The 7 - ray energy deposition is computed 
using a Monte Carlo approach inillier & Dessarlll2012l) . and for 
this study only the decay of 56 Ni is taken into account. We en- 


2 In KEPLER, the width of the mixing box is specified in mass space rather 
than velocity space. This is problematic for species like H and He that may 
be present only in the outermost layers of the progenitor star, potentially 
introducing some mixing out to the largest velocities in low mass ejecta, 
which is probably unphysical. 


force a minimum mass fraction for metals not included in the small 
network of 19 isotopes used in the pre-SN simulations, and assign 
them their so lar metallicity value . The model atoms are the same as 


those used in Dessart et al. 

2013^), with the updates for forbidden 

line transitions discussed in 

Dessart et al. ( 2014j). 


Convolving the e mergent flux computed with CMF_FLUX 
iBusche & Hillierl2005l ) over filter bandpasses and integrating over 
the full frequency range, we directly extract the photometric and 
bolometric properties of all model sequences. Because of the differ¬ 
ent physics treated in Id and CMFGEN, a few time steps are needed 
to relax the ejecta initially, in particular in the regions with a Rosse- 
land mean optical depth less than a few. Hence, we only show the 
results for expansion ages greater than m 2 d. 

The main deficiencies of the present simulations are the limi¬ 
tation to close-binary progenitors and the use of a mixing procedure 
that is simplistic. We are in the process of curing these deficiencies. 


2.2 Summary of ejecta properties 

Although model ejecta masses are lower than «5 Mq, there is 
great disparity in the ejecta composition. Model 3p65 has some 
residual hydrogen in its outermost layers, with a total mass of 
0.005 M@ and a maximum m ass fractio n of 0.14. It constitutes a 
prototype for a type lib event i Dessart et alJl20lH) . 

Models 3p65 and 6p5 have a massive He-rich shell of about 
1 Mg in which the helium mass fraction is 90%. This shell repre¬ 
sents a large fraction of the total ejecta mass, from « 70% in model 
3p65 d own to 35% in model 6p5. As we discuss in lDessart et al .1 
d2012l) . the larger that fraction the greater the likelihood of exciting 
He I lines and producing a SN lb (or lib if hydrogen is present). 
Model 5pl 1 is hydrogen deficient and exhibits a balanced mixture 
of He, C, and O at its surface. 

If all these progenitors were luminous enough to drive an 
optically-thick wind and produce an emission-line spectrum, mod¬ 
els 3p65 and 6p5 would correspond to a WN or WNh type, but 
model 5p 11 would clearly be of WC type. Model 5p 11 has lost 
its He-rich shell (that shell made of J> 90% helium). Its surface is 


© 2015 RAS, MNRAS 000101261 
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Figure 2. Bolometric luminosity light curve for models 3p65Axl. 5pl lAxl, and 6p5Axl. We use a logarithmic scale for the .r-axis in the left panel to show 
the contrast in early-time luminosity and the possible existence of a short-lived post-breakout plateau (even in the presence of mixing). The right panel shows 
the full evolution, and in particular the increase in the peak delay, light curve width, and nebular brightness for higher mass ejecta (associated in this context 
with the higher 7 -ray trapping efficiency). For each model, we overlay the decay power associated with 5 ®Ni and its daughter isotope 56 Co (dashed line). 


Table 5. Luminosity and photometric properties for our model set. For each entry, we give the rise time to maximum, the value at peak, and the corresponding 
magnitude change between the time of peak and 15 d later. The R-band maximum occurs within 1 day of the bolometric maximum. 


[mag] Model 

^rise 

[d] 

Tbol 

Max. 

[ergs -1 ] 

AMis 

[mag] 

^rise 

[d] 

Tuvoir 

Max. 

[ergs -1 ] 

AMis 

[mag] 

^rise 

[d] 

U 

Max. 

[mag] 

AMis 

[mag] 

^rise 

[d] 

B 

Max. 

[mag] 

AMis 

[mag] 

3p65Axl 

5pl lAxl 
6p5Axl 

2.631(1) 

3.584(1) 

4.212(1) 

1.648(42) 

1.590(42) 

1.495(42) 

6.157(-1) 
3.939(-l) 
2.262(-l) 

2.579(1) 

3.506(1) 

4.057(1) 

1.290(42) 

1.182(42) 

1.079(42) 

6.956(-l) 

4.234(-l) 

2.449(-l) 

2.387(1) 

3.208(1) 

3.634(1) 

-1.576(1) 

-1.524(1) 

-1.497(1) 

1.157(0) 

6.838(-l) 

2.554(-l) 

2.410(1) 

3.256(1) 

3.668(1) 

-1.662(1) 

-1.633(1) 

-1.606(1) 

1.028(0) 
5.941 (-1) 
3.623(-l) 


within the C-rich part of the CO core. This model appears more 
suitable to produce a type Ic SN. 


In this paper we focus on three models, endowed with the 
same ejecta kinetic energy at infinity of « 1.2x 10 jl erg, but differ¬ 
ent masses and composition. These models, in order of increasing 
pre-SN mass, are 3p65Axl, 5pllAxl, and 6p5Axl. The models 
have initial 56 Ni masses of 0.074 to 0.099 Mp. Thi s is som ewhat 
smaller than the mean of 0.2 Mg inferred by iDrout et ail d2011!) 
for SNe Ibc, but there is considerable scatter within the class. 
Progenitor and ejecta properties are summarised in Tables UM 
while the composition for each SN model is shown in velocity and 
Lagrangian-mass spaces in Fig. (T] 


The goal with this restricted set is to discuss the salient fea¬ 
tures of type lib, lb, and Ic rather than correlate the properties of 
their light curves and spectra with those of the progenitor and ex¬ 
plosion — this will be done in Paper II. We also address the influ¬ 
ence of mixing by comparing these models to the “x2” counterpart, 
i.e., 3p65Ax2, 5pllAx2, 6p5Ax2 (Section |5J. The description of 
the full model set is deferred to Paper II. 


3 PHOTOMETRIC EVOLUTION 
3.1 Bolometric light curve 

After shock emergence, the light curve of a SN is controlled by 
the simultaneous effects of cooling, heating, and energy transport^ 
The cooling is initially primarily caused by expansion until the 
ejecta reaches a large radius. When the ejecta turn optically thin, 
cooling is primarily by radiation. The dominant heating source is 
radioactive decay, with a small contribution from recombination at 
early times. 

The shocked envelope holds a large internal energy at shock 
breakout (about 50% of the total ejecta energy at that time). Despite 
the degradation by expansion of this shock-deposited energy, the 
SN is still significantly luminous a few days after the explosion. 
A post-breakout plateau of the order of 10 41 ergs -1 is seen in the 
larger mass model 6p5Axl, while mixing in the lower mass models 
leads to a gradual brightening by 2-3 days (Fig. 0. 


This post-breakout luminosity is function of the progenitor 
radiu s, ejecta mass, and the amount and spatial distribution of 


56 Ni (Woosley et al. 

1994; Dessart et al. 2011; Bersten et al. 201^; 

Piro & Nakadl20131: 

Nakar & Piro 20141). In this work, all oroeen- 


3 By definition, the bolometric luminosity includes contributions at all 
wavelengths. However, here we adopt the usual convention in SN model¬ 
ing, and exclude the contribution of gamma-rays which generally cannot be 
measured. 
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Table 6. Same as Table[5] but now for the V, R, and I bands. 


Model 

^rise 

[d] 

V 

Max. 

[mag] 

AM 15 

[mag] 

^rise 

[d] 

R 

Max. 

[mag] 

AM 15 

[mag] 

^rise 

[d] 

I 

Max. 

[mag] 

AM 15 

[mag] 

3p65Axl 
5pl lAxl 
6p5Axl 

2.551(1) 

3.442(1) 

3.941(1) 

-1.723(1) 

-1.715(1) 

-1.702(1) 

7.828(-l) 

4.745(-l) 

2.804(-l) 

2 . 668 ( 1 ) 

3.616(1) 

4.274(1) 

-1.734(1) 

-1.734(1) 

-1.731(1) 

6.099(-l) 
4.108(-1) 
2.656(-l) 

2.972(1) 

3.920(1) 

4.564(1) 

-1.740(1) 

-1.747(1) 

-1.754(1) 

3.855(-l) 

2.772(-l) 

2.064(-l) 


Table 7. Same as Table[5] but now for the J, H, and K bands. 


Model 

^rise 

[d] 

J 

Max. 

[mag] 

AM 15 

[mag] 

^rise 

[d] 

H 

Max. 

[mag] 

AM 15 

[mag] 

^rise 

[d] 

K 

Max. 

[mag] 

AM 15 

[mag] 

3p65Axl 
5pl lAxl 
6p5Axl 

2.751(1) 

3.786(1) 

4.655(1) 

-1.731(1) 

-1.746(1) 

-1.752(1) 

4.328(-l) 

3.279(-l) 

2.208(-l) 

3.041(1) 

3.977(1) 

4.769(1) 

-1.751(1) 

-1.772(1) 

-1.777(1) 

2.992(-l) 
2.557(-l) 
1.514(-1) 

2.862(1) 

3.668(1) 

4.610(1) 

-1.764(1) 

-1.768(1) 

-1.776(1) 

4.406(-l) 

3.500(-l) 

2.127(-1) 


itors have surface radii of < 20 Rq, which causes the faintness of 
the post-breakout phase relative to the subsequent peak. 

The SNe Ilb/Ib/Ic models systematically brighten a few days 
after explosion; this stems from energy deposited into the ejecta 
by the decay of 56 Ni exclusively (other unstable isotopes are ig¬ 
nored in this work). The three models have «0.09M@ of 56 Ni 
and radiate about the same amount of energy. They do so over 
a longer time for higher mass ejecta, and consequently reach a 
smaller maximum luminosity. For our model set, the rise time to 
bolometric maximum increases from 26.3 d for model 3p65Axl to 
41.9 d for model 6p5Axl, while the peak luminosity drops from 
1.65 to 1.49xl0 42 ergs -1 , in the same model orderQ 

The post-maximum bolometric light curve shows a more rapid 
fading for smaller mass ejecta (i.e., 0.61 mag for model 3p65Axl, 
and only 0.22 mag for the higher mass model 6p5Axl). For a higher 
ejecta mass there is more efficient trapping of stored energy (pri¬ 
marily in the form of moderate-energy “thermal” photons) within 
the ejecta, as well as more efficient trapping of the 7 -rays released 
by the decay of 56 Ni and 56 Co. The latter becomes more evident as 
time progresses. At 140 d past maximum, the bolometric light curve 
drops at 0.0196mag per day for model 3p65Axl, 0.0178 mag for 
model 5pl lAxl, and 0.0153 mag for model 6p5Axl. 


3.2 Peak bolometric luminosity and 56 Ni mass 

In our three models the peak luminosity is « 50% larger than the 
instantaneous decay rate at maximum, significantly offset from the 
prediction of Arnett’s rule (which states that the bolom etric lumi¬ 
nosity is equal to the total decay rate at that time: lArnettll982h . This 
“rule”, identified in a simplified analytical model of SN la light 
curves, is not a fundamental law that should hold exactly because, 
for example, of energy conservation. There is thus no fundamen¬ 
tal reason why it should apply accu rate ly to type Ibc SNe. In the 
delayed-detonation models of lBlondin et alj i2013h . which cover a 
range of 56 Ni masses, the offset is up to « 10%, and Arnett’s rule 


4 If we scale the peak luminosity to exactly the same 56 Ni mass of 
0.0742Mq, the peak luminosity drops from 1.65 to 1.12xl0 42 ergs -1 . 


yields a reliable estimate of the 56 Ni mass. Flowever, for the mod¬ 
els presented here, Arnett’s rule overestimates the 56 Ni mass by 
44%, 53%, and 52% for models 3p65Axl, 5pllAxl, and 6p5Axl, 
respectively. 

To estimate the 56 Ni ma ss from the bolome tric luminosity, one 
may instead use the idea of iKatz et al .1 ( I2013I) . which is robustly 
based on an energy-conservation argument. Considering heating 
from radioactive decay and cooling through expansion and radia¬ 
tion, one can solve the internal energy equation out to late times to 
find 


— tohi[to ) 


' to 


Jt 0 


where Ld ec (f) is the decay luminosity and E(t) is the total radia¬ 
tion energy trapped within the ejecta at t. At late times, as the ejecta 
becomes optically thin, there is essentially no stored radiation en¬ 
ergy so we neglect the term tE(t). To test this expression, we use 
to = 3 d. At that time, E is of the order of 10 48 erg, and is about a 
tenth of the time-integrated bolometric luminosity. Furthermore, as 
expected, the ratio 

[ dt't'Lboi{t') j f f dt't'Ldec(t') + toE(to) 

Jt 0 ' \j to 


is within a few percent of unity (a superior accuracy than achieved 
with Arnett’s rule) at 10-20 d after maximum (Fig. [3j. Neglecting 
the term to-E(fo) increases the offset by a few per cent. At early 
times, the ratio is much below unity because the decay energy is 
deposited at large optical depth, too large to influence the rate at 
which radiation escapes from the photosphere. At late times, the 
ratio drops below unity again because of 7 -ray escape — the value 
of unity would be recovered if we used the decay energy effectively 
deposited within the ejecta, or if we included the escaping 7 -ray lu¬ 
minosity in the bolometric luminosity. At 10-20 d after maximum, 
7 rays are fully trapped in typical SNe Ilb/Ib/Ic and so the relation 
holds satisfactorily without any detailed calculation at that time. 
The downside of this approach, which also affects Arnett’s rule 
(although only at maximum), is the difficulty of constraining the 
bolometric luminosity from limited multi-band photometry. 
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Figure 3. Illustration of the variation of the ratio J* dt't'L\> 0 \(t') /[ J* dt't'L^ ec (t') + toE(to)] versus time since bolometric maximum (we use to = 3 d). 
As pointed out bv lKatz et al] J201 3t> . this ratio eventually approaches unity at late times if one uses the trapped decay energy. Since the trapped fraction is not 
directly inferred, we find more convenient to check this expression using the total decay energy. In that case, the ratio is unity to within a few per cent about 
10-20 d after maximum. At later times, 7 -ray escape makes this ratio drop below unity, the more so as time passes and the ejecta thins out. Higher mass ejecta 
exhibit a longer period in which the ratio sustains its maximum value. 
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Figure 4. Multi-band light curves (left) and colour evolution for a selection of optical and near-IR bands (right) for model 3p65Axl. Redder bands tend to 
peak at later times while the colour evolution shows a complex behaviour. Notice how the minimum colour, for each plotted filter combination, occurs prior to 
/(-band (and hence bolometric) maximum. 


3.3 Multi-band light curves 

In this section, we use multi-band light curves to discuss the dis¬ 
tribution of the emergent flux with wavelength and time. Results, 
including colours, are shown in Figs. EH 6 ] while rise times, peak 
magnitudes, and post-maximum decline rates (the change in mag¬ 
nitude 15 d after the time of peak) are provided in Tables EHI] 

Generally we obtain shorter rise times in bluer optical bands, 
although an exception is that the /-band peaks later than the J 
band (likely because of the strengthening Call triplet lines; see 
Section EJ- The R band peaks within about a day of the bolomet¬ 
ric light curve maximum. As for the bolometric maximum, models 
with a larger ejecta mass peak later. Higher mass models also show 
a larger spread in the times at which each magnitude peaks. Be¬ 
cause of the comparable 56 Ni mass, the peak magnitudes are, how¬ 
ever, similar, band-by-band, for all three models. For example, we 
obtain a peak R band magnitude of —17.34, —17.34, and —17.31 


for models 3p65Axl, 5pllAxl, and 6p5Axl (a closer correspon¬ 
dence than for the bolometric luminosity maximum). 

A similar be haviour is seen in SN la models (see, e.g., 
iBlondin et al.ll20l1 ) and in both cases results from a combina¬ 
tion of effects. On the rise to maximum, the ejecta heats up and 
the emissivity at depth shifts to shorter wavelength — this is in 
part reflected by the behaviour of the electron-scattering photo¬ 
sphere (Fig.[7]>. The ejecta is thus not expected to evolve at constant 
colour. However, line blanketing causes strong energy redistribu¬ 
tion to longer wavelength. As we approach the bolometric maxi¬ 
mum, lines strengthen (see also Section |7J, dominating the emis¬ 
sivity and the opacity over continuum processes. The magnitude of 
the effect grows with time as the spectrum forms deeper into more 
metal-rich regions. 

Past maximum, the decline rates are greater in bluer optical 
bands, and the more so for lower mass models - AA /15 =1.17 mag 
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Figure 5. Same as Fig.0 but now for model 5pl lAxl. 
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Figure 6. Same as Fig. [4] but now for model 6p5Axl. The curves are similar to those shown for model 3p65Axl, but the timescales for the variations are 
longer. 


in the U band for model 3p65Axl, but only 0.25 mag for model 
6p5Axl. In the higher mass model 6p5Axl, the decline rate over 
15 d in different bands shows much less scatter than for the lighter 
model 3p65Axl, which reflects the longer diffusion time of the for¬ 
mer. 

The colour evolution for all three models is qualitatively sim¬ 
ilar, despite the differences in mass and composition (right panels 
of Figs. EJ®. During the first few days after explosion, all models 
become redder. Unfortunately, although there may be a physical ef¬ 
fect here, this may in a large part stem from the initial relaxation of 


the Id model into CMFGEN. Subsequently, because of the heat sup¬ 
ply from decay, they become bluer on the rise to maximum. From 
maximum until 20 (model 3p65Axl) to 40 d (model 6p5Axl) they 
become redder, probably because of the increasing effect of blan¬ 
keting. Beyond that time, as the ejecta turns progressively nebular, 
the colours become bluer again. At nebular times, the colours re¬ 
flect in part the distribution of forbidden-line transitions and thus 
can show great disparity between filter sets. 

In our set of models, the initial colours are systematically red, 
unlike explosions occurring in extended stars — the larger the ini- 
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Figure 7. Evolution of the ejecta properties at the electron-scattering pho¬ 
tosphere for model 3p65Axl. 



Figure 8. Comparison of the (V — R ) light curve for models 3p65Axl, 
5pllAxl, and 6p5Axl. At lOd after R-band maximum, the mean value is 
0.322 mag and the standard deviation is 0.047 mag. 


tial radius, the bluer the colo urs at early timesand the more delayed 
the phase of recombination jDessart et al JI 201 3ah . This is a bias in 
our sample of progenitors, which is limited to close-binary systems 
in w hich t he p ri mary st ar explodes with a final radius < 20 Rq. 

iDrout et aD J20 1 ill reported a small scatter of 0.06 mag about 
a mean (V — R) colour of 0.26 mag at 10 days after V -band maxi¬ 
mum for observed SNe Ilb/Ib/Ic with an inferred reddening. In our 
models, we also find that the colours early after peak are compa¬ 
rable. For (V — R) at 10 d after R-band maximum, we obtain a 
mean of 0.32 mag, with a standard deviation of 0.04 mag (Fig. H. 
The larger sample of models also shows this behaviour, with a sim¬ 
ilar mean (the standard deviation is reduced if we limit the sample 
to lower-mass ejecta; Pessart et al., in prep.). This uniformity in 
post-maximum colour can thus be used to constrain the reddening 
in SNe Ilb/Ib/Ic. 

While the bolometric correction can be estimated for our mod¬ 
els for any band, a more accurate estimate of the bolometric lumi¬ 
nosity can be obtained by using multiple passbands. Here we ex¬ 
tract the relative fraction of the flux that falls within the UV, the op- 





Figure 9. Fractional flux falling into different spectral regions for the mod¬ 
els 3p65 Ax 1, 5p 11 Ax 1. and 6p5 Ax 1 up to about 200 d after explosion. 


tical and the near-IR (Fig. [9}. Because of the cool “photospheric” 
temperature of our model (Fig. [7j- and the strong blanketing by 
lines (see below), the UV is subdominant at all times. The UV rises 
a little at nebular times, in particular in the lower mass model, but 
remains sub-dominant. The bulk of the radiation emerges in the 
optical range 80%), with the rest at longer wavelengths and pri¬ 
marily in the near-IR — little flux is emitted beyond 2-3 m. 
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X [/I m] 

Figure 10. Spectral evolution of model 3p65Axl from 0.3 to 4/i,m until 70 d after explosion (we show A 2 F\ for better visibility). Overlaid on the full synthetic 
spectrum (black thick line), we hatch the offset in flux obtained when HI and He I bound-bound transitions are omitted from the formal solution of the radiative 
transfer equation. The presence of HI and He I lines at early times, and the absence of HI at late times, makes this model a type lib. 
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A [/I m] 

Figure 11. Same as Fig. llOl but now for model 6p5Axl. HI lines are not plotted since hydrogen is absent in the ejecta. This model corresponds to a Type lb 
SN. Even for this model, with a He mass of 1.67 M©, He I lines (except for A10830 and A20580) are relatively weak and/or blended. At 70 days, and unlike 
earlier epochs, He 1 10830 A is almost entirely responsible for the P Cygni profile seen near 10830 A. 
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X [/I m] 

Figure 12. Same as Fig. EL but now for model 5pllAxl. Although helium is present in the ejecta, only the He I line at 10830 A is visible over the entire 
optical and near-IR ranges, and only at late times. This model corresponds to a Type Ic SN. At early epochs Na ID could be erroneously identified with He I 
A5876. 
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4 SPECTROSCOPIC PROPERTIES 

4.1 Spectral evolution for our model set 

The three models 3p65Axl, 5pllAxl, and 6p5Axl have the 
same explosion energy (« 1.2x 10 51 erg), a similar 56 Ni mass 
(«O.O9M 0 ), similar ejecta masses (2-5 Mg), and small progeni¬ 
tor radii (<^20Rg). However, a fundamenta l dist in ction between 
these ejecta is their composition. As shown in lDessart et all 1 201 1|) , 
even ~ 0.001 Mg of hydrogen can lead to a strong Ha line at early 
times, while even a h elium mass fraction of ^30% may produce no 
He I line jDessart et al.ll2012h . Our three models have composition 
characteristics (see Tables[3]and[4] and Fig. IT} that make them suit¬ 
able candidates for a type lib (model 3p65Axl), a type lb (model 
6p5Axl) and a type Ic (model 5pl lAxl). 

We show the spectral evolution of each model in Figs. 1 1 0TlT2l 
covering from 4 to 70 d after explosion and over the range 0.3- 
4/rm. None of our models show blue featureless spectra at early 
times, likely because the progenitor radii are small so that expan¬ 
sion cooling causes a large drain of internal energy in the outer 
ejecta layers. To avoid unnecessary speculations, we delay further 
discussion to a detailed study focused on that aspect alone. 

There is a great similarity of the spectral energy distribution 
for the three models at all times shown, with the bulk of the flux 
emerging in the optical range and strong signatures of line blan¬ 
keting by Fe II and Ti II (see below). However, closer examination 
reveals significant differences. To highlight these differences, and 
to make more obvious the classification of the model as type lib, 
lb, or Ic, we colour the flux associated with specific ions. This con¬ 
tribution is estimated by excluding the bound-bound transitions of 
a given ion in the formal solution of the radiative-transfer equation. 

Model spectra for 3p65Axl exhibit Hi lines (primarily Ha 
and H/3 in the optical, and Pa in the near-IR) at early times 
and up to the peak, while He I lines are present at all times, 
both in the optical and in the near-IR (Fig. ED- These He I 
lines may be present ea rly on without the influence of 56 Ni de¬ 
cay toessart et alj|201 ll ). but non-thermal excitation and ionisa¬ 
tion are essential to produce them after a few days jLuc 3 mu 
IPessart et al J|20 1 3 : lHachinsier et al]|2012l : pLi et all 1 201 2|) . This is 
because the temperature in the spectrum formation region is low, 
so that He I lines cannot be thermally excited. The ionisation of 
He (and of other species) in the spectrum formation region is 
low (Fig. [Hi. We find He I lines at 5875, 6678, 7068, 10830, 
11013.07, 12527.5, 12846.0, 12984.9 (and numerous other contri¬ 
butions that overlap around 1.25-1.3 /im), 15083.6 (weak), 17002.5 
(weak), 18685.4 and 19543.1 (overlap with Pa), 20581.0 (strong) 
+ 21120.2 (and other similar lines that overlap on the red edge of 
2.0581^im), 26881.2 (together with other He I lines that overlap), 
and 37025.6 A. With the presence of H I lines early on and the per¬ 
sistence of He I lines, the 3p65Axl model corresponds to a SN lib. 

Model 6p5Axl shows a very similar evolution to model 
3p65Axl but because the ejecta is hydrogen deficient, it shows no 
HI lines (Fig. Ql}. This model corresponds to a type lb SN. Inter¬ 
estingly, the two sequences for the lib and the lb models are very 
similar if we omit the differences made by HI lines. Hydrogen is 
so under-abundant in model 3p65Axl that it does not affect signifi¬ 
cantly the radiative transfer solution - hydrogen acts only as a trace 
element. 

Model 5pllAxl shows the same overall evolution in spectral 
energy distribution but for this carbon and oxygen rich ejecta the 
spectral signatures are distinct (Fig. rT2l Pespite the 0.32 Mg of 
helium and 30% surface helium mass fraction, He I lines are always 
absent in the optical. The He I line at 10830 A is present at all times, 
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Figure 13. Evolution of the ionisation fraction of He in model 3p65Ax 1 
(top) and of O in model 5pl 1 Axl (bottom). Specifically, the quantity plotted 
is y~\ iNj/Ntot, where Ay is the species number density when /-times 
ionised and Ntot is the species total number density. Symbols denote the 
location of the photosphere (using the Rosseland-mean opacity). The times 
shown are 4, 6, S, 10, 15. 20. 30, 40, and 50d after explosion for the top 
panel. The times shown are 5, 10, 15,20, 30,40, 50, and 70 d after explosion 
for the bottom panel. 


but it is badly blended at early times - absorption that could be 
erroneously associated with He I alone is also due to C I and Mg II. 
Paradoxically He I A10830 is strongest at late times (see below). 
This model is significantly mixed (just like the other two models) 
and so the lack of He I lines res ults fr om t he low helium abundance 
and the large CO core mass JPessart et al |20 1 3l . 

The buffer of mass lying in between the 56 Ni rich layers and 
those rich in helium is much larger in model 5pllAxl, prevent¬ 
ing efficient non-thermal excitation at early times. Model 5pl 1 Axl 
exhibits an extended oxygen-rich shell, whose outer part is rich in 
carbon and poor in helium (mass fraction of 10%). In contrast, in 
models 3p65Axl and 6p5Axl, the oxygen shell is less massive and 
resides underneath a massive helium-rich shell (>90% mass frac¬ 
tion helium). This difference in core structure stems from the dif¬ 
ference in main sequence mass. Models 3p65 and 6p5 stem from 
stars with initial masses 16 and 25 Mg, which make a relatively 
small CO core compared to model 5p 11 which stems from a star 
with an initial mass of 60 Mg that eventually formed a massive CO 
core. All else being the same ( 56 Ni mass, mixing magnitude, ex- 
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Figure 14. Maximum light optical spectra for models 3p65Axl (SN lib), 6p5Bxl (SN lb), and 5pl 1 Axl (SN Ic). We illustrate the contributions from lines of 
different ions by colouring the space that separates the full spectrum (dark black curve) from that obtained by neglecting a given ion from the formal solution 
of the radiative-transfer equation. 
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Figure 15. Same as Fig. 1 141 but now for the near-IR range. For better visibility, we show the quantity A 2 F\. 
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plosion energy), this difference in core structure causes the lack of 
He I signatures in model 5pl lAxl. In this SN Ic model, the feature 
at 5900 A is solely due to Na ID (strong at early and late times but 
weak around maximum). 

Our result disagrees somewhat with lHachinger et al .1 J2012h 
who find in their study that as little as 0.14 Mq of He may suffice 
for the production of He I lines. There is probably no clear thresh¬ 
old because the excitation of He I lines will depend not just on He 
abundance, but also on 56 Ni mass, mixing, the size of the CO core, 
and the C and O mass fraction in the outer parts of the progenitor. 

Iron group elements have a mass fraction equal to the progen¬ 
itor metallicity value throughout the ejecta except in the innermost 
layers (with the exception of j6 Ni and its daughter nuclei). Here, 
iron-group elements cause blanketing in the optical (Fe II, Ti II) and 
in the UV (Ti II). Numerous lines from Co II produce a broad emis¬ 
sion feature at 1.6/im. Because of steady burning within the pro¬ 
genitor helium core, intermediate mass elements have abundances 
that are over- or under-abundant compared to solar, with variations 
between models and also with depth for each model. 

In the type Ic model, most of the features between 9000 A and 
« 1/rm are due to C I, while these features are very weak in the 
Type lib model 3p65Axl. In the type lb model 6p5Axl, these CI 
lines are visible around maximum and later because of the larger 
CO-rich shell in this more massive progenitor. Hence, CI lines offer 
a possible diagnostic if the identification of He I lines is ambiguous 
- a strongly He-rich photosphere would most likely exhibit weak 
or no CI lines in the 9000 A and « 1 /.tm region. 

Oxygen shows only one signature in the optical with 
01 7777 A. This feature is strong at all times in model 5pllAxl 
but its presence in the type lib model 3p65Axl likely results from 
mixing (the helium-rich shell is strongly deficient in oxygen in the 
pre-SN model). Improving our current algorithm for mixing is work 
in progress — while 56 Ni may b e mixed far out in the ejecta, inter- 
mediate mass elements may not jWongwathanarat et all2014h . 

Na ID is conspicuous only in the type Ic model, and could be 
confused with He 15876 A. Neon lines are absent at all times and in 
all models. Magnesium lines stem primarily from Mg II. In the opti¬ 
cal, we predict (usually blended) lines at 7790, 8213-8234, 9218— 
9244, 10910, 21368-21432, 24042-24124 A. In model 3p65Axl, 
these lines are probably present even at early times because we mix 
magnesium out in velocity space (O, Ne, Mg are originally under- 
abundant in the He-rich shell). 

Lines from Call (H&K lines and the triplet at 8500A) are 
strong at all times, and of comparable strength in the three mod¬ 
els. Being closely tied to the ground state (directly or through 
Call H&K for the triplet), these lines are strong even at solar abun¬ 
dance so they provide little constraint to help distinguish different 
ejecta. 

All these signatures are illustrated at the time of bolometric 
maximum for the optical range in Fig.[l4]and for the near-IR range 
in Fig. ED These figures are merely for illustration, not quantitative 
estimates, since excluding a given ion in the computation of the 
spectrum breaks the consistency of the calculation. 


5 NEAR-IR He I LINES 

The identification of He I lines in SN Ib/c spectra is a recurring 
issue in the SN community. The problem is overemphasised since 
helium is likely present in all core-collapse SNe, simply in different 
abundances. For example, optical He I lines are seen in Type II SN 
spectra only prior to hydrogen recombination. During the rest of 


the photospheric phase, optical He I lines are absent a lthough ab out 
30% o f the total hydrogen-envelope mass is helium l IDessart et al.l 
l2013bl) . In that case, there is little debate about the presence of 
helium because hydrogen is present, although the difficulty in pro¬ 
ducing He I lines in such cool atmospheres is already evident. 

In hydrogen poor ejecta, the same issue returns, but here the 
problem is whether the SN atmosphere lies within the He-rich shell 
(>90% helium mass fraction) or within the CO core (typically 
10% helium mass fraction). Hence, the production of He I lines 
depends on the ejecta composition. Furthermore, the need for non- 
thermal excitation/ionisation represents a major hurdle. Any ejecta 
with a massive CO core will probably never show strong He I lines 
because of the large mass buffer between the 56 Ni rich regions and 
the outer He-rich regions (whose presence is not even guaranteed 
because of mass loss). This suggests that SNe Ilb/Ib progenitors 
preferentially arise from lower mass progenitors (which have small 
CO cores), where ev en mo derate mixing can foster the production 
of He I lines dDessart et al] 2012 l) . 

In this section, we document the properties of two spectral 
regions notorious for their possible He I signatures, at 10830 A 
and 20581 A (Figs. [ 1611181 1 . He I 10830 A is the strongest line, and 
has the pot ential to confirm or deny the presence of He in SN Ic. 
iFilippenko et al.lll995l ) used the presence of strong absorption near 
10830 A to ar gue for the pre sence of He in the Type Ic SN 19941. 
Later work bv IMillard et al.l d 19991 ) found, using the parameterised 
synthesis code SYNOW, that the feature could not be fitted by He I 
alone. They sug geste d that it was a possible blend with CI. Later 
work by Sauer_etak[d200^ confirmed this sugges tion. Other stud¬ 
ies (e.g., Meikle et a l.ll996i : lHachinger et al.l2012h show that Mg II 
is also a possible contributor. 

In type Ilb/Ib models 3p65Axl and 6p5Axl, the 10830 A re¬ 
gion is dominated by the He I line, with a modest contribution of 
Mg II10914 A (multiple transitions), and 0111300 A in the red 
wing (multiple transitions). At late times, numerous lines from 
Call, Fe II, and CoII appear around Ifim and affect the blue wing 
of the He I line. In model 5pllAxl, similar features are present, 
but with the addition of a strong overlapping component from nu¬ 
merous Cl transitions around 10700 A. Around bolometric maxi¬ 
mum, the He I contribution to the 10900 A feature is weak relative 
to that of CI and Mg II. However, it is still important. To remove the 
feature seen near 10900 A we would need to simultaneously omit 
He I, C I, and Mg II lines from the spectrum calculation. The rela¬ 
tive strength of these different contributions varies with time, and 
He I is a significant component at all times only in the type Ilb/Ib 
models. 

In the 20581 A region, the dichotomy between SNe Ilb/Ib and 
SNe Ic is more evident. Indeed, He 120581 A is present in our type 
Ilb/Ib models but completely absent in the type Ic model. Because 
there are no other lines at that location, the lack of a feature th ere is 
clear evidence that there is no He I line. lHachinger et all d2012l ) also 
highlight the importance of the He 120581 A line in assessing the 
He abundance of SN Ibc ejecta. Note that in our SNe Ilb/Ib models, 
there are He I lines around 1.9 /.tm (e.g., 19543.1 A), which spread 
the influence of He I far to the blue of the 2.0581 gm feature. 

An i nteresti n g fea ture of our simulations, also observed (see, 
e.g.. lErgonet"ai]|2014l for SN 201 ldh), is the migration of the lo¬ 
cation of maximum absorption (or more generally the broadening) 
of the He I lines at 10830 and 20581 A. This is the spectral counter¬ 
part of the rapid photometric decline at late times. Indeed, before 7 
rays start to escape, they first deposit their energy non locally within 
the ejecta, increasingly influencing the outer layers after bolometric 
maximum. When this effect arises, non-thermal excitation and ion- 
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Figure 16. Illustration of the 10830 and 20581 A regions for model 3p65Axl (type lib). We show the contributions from He I, C I, O I, Mg II lines to highlight 
the fractional importance of each species/ion. In these lib spectra, high velocity absoiption (around —10,000 km s — - 1 ) associated with both 10830 and 20581 
provides unambiguous evidence for the presence of He. However, even in this case Mg II has a significant influence on the emission component of 10830 A. 




Figure 17. Same as Fig. ED but now for model 6p5Axl (type lb). He 120581 A is present at all times, and hence provides a very useful diagnostic of the 
presence of helium. 




Figure 18. Same as Fig. DU but now for model 5pl 1 Axl (type Ic). Notice the complicated blend associated with the P Cygni profile at 10830 A while there is 
a firm non-detection of He I at 20581 A. At late times the 10800 A feature shows high velocity absoiption (beyond —10,000 km s —1 ) that is only associated 
with He 1 10830 A. 
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Figure 19. Ejecta composition profiles for models 3p65A (top) and 5pllA 
(bottom) with moderate (xl) and strong (x2) mixing. [See Sections [2]and[6] 
for discussion.] 

isation processes are boosted in the outer ejecta layers (which are 
poor in 56 Ni initially), increasing the optical depth in these strong 
He I lines. In type Ic model 5pllAxl, only the He I contribution 
to the 10900 A feature shows this behaviour, which may serve as a 
tracer of helium in the ejecta. 


6 INFLUENCE OF MIXING 

In lDessart et al.1 ( l2012h . we performed controlled experiments on 
the effect of mixing in He-rich low-mass ejecta models. This study 
emphasised how the distribution of 56 Ni is critical for the produc¬ 
tion of He I lines. Indeed, even in ejecta abundant in helium, weak 
mixing would lead to a type Ic SN. The existence of SNe Ilb/Ib is 
a robust indication that extensive (macroscopic) mixing takes place 
in these explosions, and by extension, is a likely feature of all core¬ 
collapse SNe0 

The production of He I lines and a type lb SN requires both 
the presence of a large helium mass and a strong mixing of 56 Ni 

5 We note that, theoretically. He I lines may be present at the earliest times 
without the influence of non-thermal effects, but their presence in this case 
is very short lived and limited to progenitors with a He-rich shell (for a 
discussion, see[Dessart et al.l 201 ill . 


Figure 20. Top: Illustration of the effect of a variation in the magnitude of 
mixing (xl versus x 2 , i.e., moderate versus strong mixing) on the bolomet- 
ric luminosity light curve for our three ejecta models. Bottom: Same as top, 
but now showing the bolometric light curves normalised to their maximum 
value and shifted so that the origin is the time of bolometric maximum. 

into the He-rich layers. This implies that any progenitor star char¬ 
acterised by a large CO core would tend to form a SN Ic because 
the He-rich regions are too remote from the inner ejecta where 56 Ni 
is produced. It is likely that mixing is a ubiquitous feature of core¬ 
collapse SNe. The details and degree of mixing as a function of 
progenitor and explosion properties are, however, uncertain. 

In the present work, all models have therefore been mixed. 
From our set of models, we indeed produce type Ilb/Ib and Ic 
SNe for the reasons described above. But to quantify the uncertain¬ 
ties associated with mixing, we run strongly-mixed counterparts to 
models 3p65Axl, 5pllAxl, and 6p5Axl — these strongly mixed 
models are given the suffix x2. Our mixing procedure at present ap¬ 
plies to all species and is thus indiscriminate (Fig. Q9>. In contrast, 
models without mixing would show shells with nearly 100% 56 Ni 
mass fraction in the innermost layers of the ejecta and non e above, 
as in the (unmixed) models studied i mPessart et al.1 j2012f) . 

As is well known, s tronger mixing cau ses a higher luminos ity 
(e.g jBersten et al.ll2012l) and bluer colours l lDessart et alj|2012t) at 
early times. The larger amount of 56 Ni at larger velocities causes 
enhanced escape as time progresses, causing a lower luminosity af¬ 
ter maximum. In all three models, mixing causes the maximum to 
shift to earlier times, with the largest shift occurring for the high- 
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Figure 22. Illustration of the impact of mixing on the spectral evolution of our three models. Each pair corresponds to ejecta with a moderate (xl) and a 
strong (x2) mixing. Labels provide the time since explosion as well as the magnitude offset between the two models at each date. Besides the impact on the 
bolometric luminosity and colour, stronger mixing boosts non-thermal processes and leads to broader lines (at fixed E^ n and M e ). The effect is therefore 
more evident on He I lines. 



-40-20 0 20 40 60 80 100 -40-20 0 20 40 60 80 100 


Days since 7?-band Maximum 

Figure 21. Illustration of the impact of mixing on the colour evolution of 
model 5pl 1 A. Time is given in days since R-band maximum. Model with 
suffix xl has weaker mixing than model with suffix x2 (see Fig. Q9] for the 
composition profiles). 

est mass ejecta (Fig. [20t . For 3p65Axl/x2 the shift is 1.4 d, for 
5pllAxl/x2 it is 5.0d, and for 6p5Axl/x2 it is 6.7d. For model 
6p5Ax2 the maximum bolometric luminosity is about 10% larger 
compared to the model with the lower mixing; for the other models 
the change is even smaller. When we use the time of bolometric 
maximum as the time origin, the curves for the moderately mixed 
and strongly mixed models show a remarkable similarity for the 
two models with the lowest ejecta mass - the bolometric light curve 
for 6p5Ax2 is somewhat broader (by about a week at 50% maxi¬ 
mum) than the light curve for 6p5Axl (Fig. l20t . This morpholog¬ 


ical similarity for different mixing means that the explosion date 
could be in error by as much as a week if only the bolometric light 
curve was used. 

In our procedure, mixing enhances the mean atomic weight of 
the outer ejecta layers. While this should strengthen line blanket¬ 
ing, we observe bluer colours in strongly mixed models because of 
the extra heating from 56 Ni decay (Fig. ED- After maximum, the 
colour offset due to different mixing is much weaker. For (V — R) 
at 10 d past i?-band maximum, the offset is negligible. 

Synthetic spectra corroborate these properties but also il- 
l ustrate how stron ger mixing leads to broader lines (see also 
iDessart et alJl20l3) . even at fixed ejecta mass and kinetic energy 
(Fig. 1221) . The effect is particularly visible on Fie I lines and di¬ 
rectly stems from the modulation of the magnitude of non-thermal 
processes — SN Ic spectral line profiles are less sensitive to mixing 
in that respect. Uncertainties in the mixing magnitude in SNe lb in¬ 
troduces an uncertainty in the inferences on the explosion energy, 
ejecta mass etc. 

The impact of mixing is greater in higher mass ejecta models 
like 5pl lAxl and 6p5Axl because mixing is the only way the outer 
ejecta can be influenced by non-thermal processes. In lower mass 
ejecta, like model 3p65Axl, strong mixing is less critical because 
the outer ejecta (basically the He-rich shell) is close to the 56 Ni rich 
layer - it is only separated by a thin O-rich layer. 


7 THE AMBIGUOUS NOTION OF A PHOTOSPHERE 

For stellar work, the photosphere is usually defined to be the lo¬ 
cation at which the Rosseland-mean optical depth is 2/3, a def¬ 
inition which follows from the Eddington-Barbier relation for a 
plane-parallel atmosphere (e.g. jMihalaa 19781) . For extended atmo¬ 
spheres a more practical relation might be the atmospheric location 
where the probability of photon escape is 50%, or alternatively, that 
location beyond which 50% of the observed flux originates. Unfor¬ 
tunately, because the opacity has a strong wavelength-dependence. 
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and because of the rapidly expanding ejecta, the photosphere is not 
uniquely defined, and even at a single wavelength photons will es¬ 
cape over a region which is spatially extended. The strong effect 
of line blanketing in the UV increases the effective photospheric 
radius for bluer photons, while redder photons in the optical, de¬ 
couple from the gas at a radius comparable to the Rosseland-mean 
or flux-mean photosphere (Fig. 1231) . Still, even within the optical, 
this decoupling radius can vary by a factor of 2-3. 

In type II-Plateau SNe, the continuum opacity dominates at 
early times throughout most of the optical, and even at late times, 
there are fairly clean “continuum” regions devoid of line opacity 
and emi ssivity, for example, between Ha and the Call triplet at 
8500 A dDessart & Hillierll2005h . In type Ilb/Ib/Ic, the large abun¬ 
dance of He (and its low ionisation in the photospheric layers) and 
intermediate mass elements reduces the importance of continuum 
opacity sources and favours the role of lines (Fig. 1241 . Electron 
scattering is sub-dominant and may represent only a third of the 
total Rosseland-mean opacity, even at bolometric maximum. 

As time progresses and the ejecta turns nebular, the influence 
of lines grows. In our simulations, this causes the steep decline of 
the continuum flux below the representative flux level of the model, 
which only shows a quasi-continuum (Fig. [25}. The emergent spec¬ 
trum even at the peak is a collection of overlapping lines, causing 
simultaneous absorption, scattering, and emission. P-Cygni pro¬ 
file formation is altered from the general conception since it is not 
continuum photons that interact with the line but instead the back¬ 
ground radiation arising from overlapping lines. 


8 INFERENCE OF T HE EJECTA EXPANSION RATE 

One important characteristics of SN ejecta is their kinetic energy. 
It represents the left-over explosion energy after unbinding the pro¬ 
genitor envelope. To infer the ejecta kinetic energy, we need to infer 
its expansion rate. It is customary in the community to do this by 
constraining the location of the photosphere. 

However, it is clear from the preceding section that photons 
subject to continuum opacity will decouple much deeper in the 
ejecta than photons overlapping with a strong line, such as Ha in a 
SN lib at early time. To different lines correspond a whole range of 
intermediate heights and legitimate but distinct photospheres. 

Determining the location of these photospheres is nonetheless 
an interesting exercise. Using the location of maximum absorption, 
one can constrain the chemical stratification of an ejecta, bearing 
in mind that a line profile is also subject to multiple processes and 
uncertainties (e.g., ionisation, details of atomic physics, blending). 
In Fig. [26] we show a spectral sequence in the optical for model 
3p65Axl. Using synthetic spectra for individual species, we can 
unambiguously determine the location of maximum absorption in 
lines of sufficient strength and not overlapping with other lines of 
the same ion. We show the results (one red mark for each line and 
each epoch) for Ha, He I lines at 5875.7 A, 6678.1 A, 7065.2 A, 
and 10830.2 A, 017773.0 A, and the Call triplet at 8569.0 A. We 
omit the lines from Fe II or Ti II in the blue part of the optical since 
their contributions overlap — we see broad absorptions with inter¬ 
nal dips which prevent the unambiguous identification of a specific 
line. 

The guidance of the model helps for the velocity measurement 
on the weak optical He I lines at early times, or the vanishingly 
small line strength of Ha at and beyond bolometric maximum. It 
helps also with ambiguous dips. The depression at 8000 A in the 
spectrum at 60.6 d is not caused by Ca II absorption at high velocity 
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Figure 23. Wavelength variation of the photospheric radius R(r = 2/3), 
normalised to the radius of the electron-scattering photosphere -R(r eS ec = 
2/3), when all opacity sources are included, and when we use the 
Rosseland-mean opacity kr, oss or the flux-mean opacity ttfl ux . From top 
to bottom, we show the results for model 3p65Axl at 10. 21. and 41 d after 
explosion. The distinct locations traced by each curve emphasise the ambi¬ 
guity of the notion of a photosphere in these Type Ilb/Ib/Ic SNe. 
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Figure 24. Variation of the mass absorption coefficient associated with the 
Rosseland-mean opacity kr oss (solid) and the electron-scattering opacity 
K e sec (dashed) in the ejecta of models 3p65Axl (SN lib), 6p5Axl (SN 
lb), and 5pllAxl (SN Ic) at «4d after explosion (top) and at bolometric 
maximum (bottom). The same ordinate range is used in each panel. All 
opacities show a large variation with ejecta location, reflecting variations 
both in composition and ionisation (see also Fig EJ. The range in mass 
absorption coefficient is significantly greater than for Ty pe la SNe, which 
posse ss a more uniform composition (see, e.g.. Fig. 12 in lDessart & Hillien 


but corresponds to the blue edge of its emission. It also allows to es¬ 
tablish with certainty that the progressive migration of the 1.05 gm 
absorption to larger velocity is associated with the broadening of 
the He 1 10830 A line, following the growth of the 7 -ray mean free 
path and the non-thermal effects on He I in the outer ejecta. 

Repeating the process for the other two models, we illustrate 
the trajectory of the velocity at maximum absorption for repre¬ 
sentative lines in Fig. [27] (we use the gf -weighted mean wave¬ 
length for multiplets). In the type lib model, v a bs(Ha) converges to 
« 13000 km s~ , which indeed corresponds to the base of the H- 
rich shell in the ejecta (Tabled. Some sudden variations can occur 
when line profiles have double dips of comparable strength or when 
the absorption is an extended trough. The jumps at early times may 
also arise when we switch on the treatment of non-thermal pro¬ 
cesses. However, overall, the evolution of all lines is to exhibit a 
maximum absorption that recedes in velocity space as we approach 
the time of bolometric maximum. 


After maximum, non-local energy deposition and non-thermal 
processes alter the line opacity in the outer ejecta and lead to 
the broadening of the absorption troughs. The effect is strong in 
He 1 10830 A, whose maximum absorption occurs at 12000 km s ” 1 
at late times. The chemical stratification matters here since He is 
systematically more abundant in the outer ejecta of these models — 
intermediate mass elements are less likely to exhibit this behaviour 
because their abundance peaks in the inner ejecta. 

While the photosphere, and its associated velocity, are poorly 
defined, the representative ejecta velocity v m = y/2Ekin/ M e is 
a fundamental characteristic of the ejecta and an important pa¬ 
rameter of analytical models for SN light curves. In our three 
simulations, we find that around bolometric maximum, this quan¬ 
tity is matched to within « 10 % by t; a bs(He 1 5875 A) in model 
3p65Axl (Vm =7476 km s -1 ), or v a bs(017773A) in models 
6p5Axl (v-m =5047 km s -1 ) and 5pllAxl ( v m = 5959 km s -1 ). 
This correspondence is useful to constrain the ejecta kinetic energy 
and mass and drops the additional step, with its associated ambigu¬ 
ities, to locate “the” photosphere. In Paper II, we will study how the 
quantity v m correlates with various lines in the full grid of models. 
This exploration also needs to be extended to progenitors with large 
surface radii. 


9 CORRECTION FACTORS 

Comparing a SN spectrum to a blackbody can give some infor¬ 
mation on the colour temperature. In type II SNe, the comparison 
is meaningful at early post-explosion times because there are ex¬ 
tended continuum regions where the flux resembles a blackbody, 
merely “diluted” through the effect of scattering in the atmo sphere 
dMihalaslll978l:lEastman et~aI1ll996l:lD~e.ssart & HiHiel2005l) . 

One must apply a correction to this blackbody spectrum in or¬ 
der to obtain the total luminosity and to infer the distance to a type 
II SN. Unfortunately, in practice, one is led not so much to correct 
for the effec t of scattering but inste ad for the effect of lines (see 
discussion in lDessart & Hillierl2005i) . With type Hb/lb/lc SNe, it is 
also possible to produce such correction factors, but the correction 
arises primarily from the effect of line emission/absorption. The 
second problem, as discussed above, is that there is no unambigu¬ 
ous photosphere, since it depends much on the opacity considered. 

None theless, we have repeated the calculations of 
iDessart & Hillierl 12005b (this study was based on steady state 
models; correction factors for time-dependent simulations of type 
II SNe was discussed in IDessart & Hillien 120080 on the present 
set of models. Wejollow the procedure described in Section 2 of 
IDessart & Hillierl l i2005 ). and specifically compute the correction 
factors through a minimisation of their Equation 9. Since electron 
scattering is not a good representation of the total opacity, we 
employ the Rosseland mean instead and focus on times when the 
ejecta Rosseland-mean optical depth is > 10 (i.e., the photospheric 
phase). Velocities associated with the photosphere, and how to 
infer them, can be gleaned from Fig. 1271 

Correction factors show a complex behaviour (Fig.l28t. Unlike 
type II SNe during the photospheric phase, which show a mono¬ 
tonic evoluti on towards cooler tempe rature at fixed photospheric 
composition iDessart & Hillierl 1201 lh . type Ilb/Ib/Ic SNe exhibit 
a non-monotonic colour evolution (Fig. EM} with a photospheric 
composition that becomes progressively metal dominated. This is 
the reason for the bifurcation in the distribution of correction fac¬ 
tors (compare the top row panels with the bottom row counterparts). 

When we limit the time span to times around bolometric max- 
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Figure 25. Evolution of the total flux (black) and the continuum flux (blue) in the optical and near-IR for the models 3p65Axl, 5pl lAxl, and 6p5Axl during 
the photospheric phase. For better visibility of the near-IR, we show the quantity A 2 F\. For the continuum flux calculation, only bound-free and free-free 
processes are included in the formal solution of the radiative-transfer equation. Note the growing separation between the total flux and the continuum flux as 
time progresses, caused by the increasing role played by lines over pure continuum processes. At many epochs no pure continuum region can be seen. In the 
optical, lines act through both emission and absorption, and in a very tangled way, while at early times in the UV lines cause a flux deficit through a blanketing 
effect. 


imum, the distribution of points shows more uniformity (bottom 
row of Fig. 128b. For convenience we have fitted the correction fac¬ 
tors using a second-order polynomial in the inferred colour tem¬ 
perature. The polynomial coefficients are provided in Table [ 8 ] We 
also note that the correction factors are systematically offset be¬ 
tween models, likely arising from the distinct composition in the 
atmospheres of each model. 

These correction factors have to be used with circumspection. 
In Fig[29] we show the blackbody fits to the different optical bands 
for model 3p65 Axl around bolometric maximum. The present con¬ 
text is far from the original notion of flux dilution in scattering at¬ 
mospheres. The correction is instead used to recover from the ques¬ 
tionable approximation of the model flux as a blackbody. 


10 CONCLUSIONS 

This paper is the first of a series describing a set of radiative transfer 
simulations for SNe lib, lb, and Ic from binary star progenitors. 
In this paper we focused our attention on three models that were 
selected for their distinct composition. Model 3p65 A has a sizeable 
He-rich shell and some residual H in the outermost ejecta layers. It 
produces HI lines early on, and He I lines at all times, making it 
a Type lib SN. Although more massive, the ejecta associated with 
model 6p5A is hydrogen deficient and has 1.65 Mq of He (« 35% 
by mass). It produces He I lines and is of Type lb. Model 5pllA 
stems from a typical WC star composition, with a mixture of He, 
C, and O in the outermost ejecta. Even with strong mixing, it shows 
no He I line in the optical and is thus of Type Ic. 

Despite these differences in composition, which produce dif¬ 
ferent SN types, the qualitative properties of the bolometric light 
curves, multi-band light curves, and colours are comparable for all 
models. Having the same ejecta kinetic energy, the differences stem 


primarily from the different ejecta masses. As expected, the higher 
the ejecta mass, the longer the rise time, the broader the light curve, 
and the slower the post-maximum decline rate. These properties 
stem from the greater trapping efficiency of low-energy photons 
(which in number make the bulk of the ejecta internal energy) and 
of high-energy photons (which originate from radioactive decay). 
About 80% of the total flux falls within the optical range, with the 
rest in the near-IR — little flux comes out in the UV except at very 
late times. At early times this holds because our progenitor mod¬ 
els have small radii — it is not a fundamental property of all SNe 
Ilb/Ib/Ic. 

Quantitatively, Arnett’s rule does not produce a very accurate 
prediction of the 56 Ni mass - applied to our models it leads to 
an overest imate of the 56 N i mass by about 50%. In contrast, the 
method of lKatz et alj J2013I) works very accurately with our mod¬ 
els. Using an expression based on a solution of the energy equation, 
we can recover the original 56 Ni to within 5% by integrating the 
bolometric luminosity over time starting at 3 d after explosion (this 
offset can be reduced to <J2% if we take into account the trapped ra¬ 
diation energy at 3 d). The challenge with this method is the need to 
monitor the multi-band brightness from very early times (the later 
we start, the more inaccurate is the method), while a better accu¬ 
racy requires the knowledge of the trapped radiation energy at the 
start of the integration. 

As observed, we find that our models exhibit an earlier light- 
curve peak for bluer bands. Despite the uncertainties in mixing and 
ejecta mass/composition, the early post-maximum colours are com¬ 
parable. The (V — R) colour shows a good uniformity with a value 
of 0.32 mag (standard deviation of 0.04 mag). The effect of mixing 
is to enhance the brightness and make the SN bluer at early times. 
This leads to a flatter slope of the bolometric light curve as it ap- 
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Figure 26. Spectral sequence for model 3p65Axl during the photospheric 
phase. We overlay the location at maximum absorption in Ha, in the He I 
lines at 5875.7A, 6678.1 A, 7065.2 A, and 10830.2 A, in 017773.0 A, and 
in the Ca II triplet at 8569.0 A. To better reveal the evolving velocity shifts 
for each transition, we also draw a vertical line (dashed and coloured) 
through the maximum absorption closest to the rest wavelength. With most 
features, the location of maximum absorption shifts to smaller velocities 
with time. However for He 110830 A this is only true initially - at late times 
the maximum absorption shifts to larger velocities due to the influence 7 - 
rays and non-thermal effects in the outer ejecta. 


proaches maximum. The effect of mixing on the multi-band light 
curves is only of minor importance after maximum. 

We emphasise the spectral differences between our three mod¬ 
els. We find that Hn is eventually a good tracer of the velocity 
of the outer H-rich shell in a SN lib. We also reproduce the non¬ 
monotonic behaviour of He I lines with time, which tend to become 
narrower from early time to bolometric maximum, before becom¬ 
ing broader again after maximum. This property is a vivid illus¬ 
tration of the combined effects of non-local (decay) energy depo¬ 
sition and non-thermal processes. This effect is primarily seen in 
He I lines, which are the most sensitive to non-thermal processes in 
these SNe. 

While He is present in all our models, it can be difficult to 
identify in spectra. Due to blending, the 1.1 gm region must be 
used with caution when trying to determine whether the ejecta con¬ 
tains He. In our type Ic model lines due to Cl and Mg II are very 
important contributors to the P Cygni profile seen near 1.09 /im, 





Days since Bolometric Maximum 


Figure 27. Evolution of the velocity at maximum absorption in represen¬ 
tative lines for the three models 3p65Axl (SN lib), 6p5Axl (SN lb), and 
5pl lAxl (SN Ic). We use the ^/-weighted mean wavelength for multiplets. 
We also add the photospheric velocity (for which we employ the Rosseland- 
mean opacity). The dashed line corresponds to the representative ejecta ve¬ 
locity given by yj2E^ n /M e . 
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while in the type lib model Mg II is an important contributor. Other 
lines, such as He I 2.0581 fi m, can be used to confirm the presence 
of He I in spectra, although the absence of this line does not pre¬ 
clude the presence of A1O83O0 Paradoxically, in both our lb and 
Ic models, He 1 10830 A is strongest at late times (« 70 d). Another 
method of distinguishing the type of progenitor is to use CI lines. 
We obtain weak or no C I lines in the 9000-11000 A region prior to 
maximum in our SNe Ilb/Ib models, which have a sizeable He-rich 
shell, while the SN Ic model, with its outer layers of He, C, and O, 
shows strong C I lines in this region. 

Because a mass buffer rich in carbon and oxygen separates 
the helium-rich shell from the original site of 56 Ni, the production 
of He I lines in our models requires mixing. Without mixing, none 
of our models ar e of typ e lb (or lib), but would instead be of type 
Ic (or lie) — see lDessart et alj d2012l ) for discussion. Furthermore, 
varying the mixing strength from moderate to strong (models xl 
and x2) did not alter the SN type of these models. As done for ob¬ 
servations, we discuss the SN type of each of our models based on 
its spectroscopic properties. However, it is clear from the character¬ 
istic width of their light curves that models 5pl 1A and 6p5A have a 
too low Ekin/Me to match SN Ibc observations (see, for example, 
iDrout etal.ll201ll) . 

As for SN la, a true continuum is only seen at very early 
times. After the spectrum formation zone recedes into the metal 
rich layers of the ejecta, metal line opacity and emissivity com¬ 
pletely dwarf continuum processes. Line overlap then conspires to 
produce a quasi continuum. Consequently, the mass absorption co¬ 
efficient is greater than that due to electron scattering, and varies 
with ionisation — it is not a constant in time or in depth. Further¬ 
more, the ubiquity of lines in the spectrum formation region makes 
the notion of a well-defined photosphere obsolete. 

Despite the complexities of spectrum formation, useful in¬ 
ferences about the ejecta can still be made. For example, the lo¬ 
cation of maximum absorption in various lines (e.g., Her or He I 
lines) can serve as a diagnostic of chemical stratification. We also 
discussed how the mean ejecta expansion rate, defined by v m = 
277kin/(He correlates with the P Cygni profile absorption veloc¬ 
ity for various lines around the time of maximum. With a measure 
of both the mean expansion rate, and the ejecta mass, the energy of 
the SN explosion may be inferred. 

Finally, we computed the correction factors to invoke when 
approximating the synthetic spectra with a blackbody. Because SN 
colours are affected by line opacity and emissivity, the correction 
depends on the selected band and on the adequacy of the model for 
a particular observation. We find that our models compare well to 
observations around the time of maximum (see Paper III) so the use 
of these correction factors should be limited to such epochs. 
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Table 8. Coefficients of the second-order polynomial used to fit the distribution of correction factors for different optical-band combinations and for models 
3p65Axl, 5pllAxl, and 6p5Axl over the time span [fmax/5, l.lxtmax ]. The polynomial is of the form £(T) = y^\_ n aj{\/T^Y , where T 4 is the 
temperature in units of 10 4 K. 
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Figure 28. Top row: Distribution of correction factors for models 3p65Axl, 5pl lAxl, and 6p5Axl for optical band combinations S = ( B , V) (left), ( B , V, I) 
(middle), (VI) (right), and based on photospheric-phase epochs. Each set S must contain at least two bands to constrain the colour temperature. Bottom row: 
Same as top, but now limiting the epochs to the range [fmax/5,1.1 x£ max ], where £ max is the time of bolometric maximum. Also plotted is a fit to the 
distribution of correction factors using a second-order polynomial. The coefficients of the fitted polynomial are given in Table [8] The correction factors and 
color temperatures shown here are tabulated in the appendix. 





Figure 29. Comparison, near the time of bolometric maximum, of the total and continuum fluxes for model 3p65Axl with the blackbody flux computed with 
the associated correction factor and colour temperature Ts (see Table I All in the appendix). From left to right, we show the results for the sets S = (B, V), 
(B, V, I), and (V, I). 
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APPENDIX A: CORRECTION FACTORS FOR THE 
THREE REFERENCE MODELS 

In this section, we tabulate the correction factors for each of the 
three models during the photospheric phase. These SNe Ilb/Ib/Ic 
models correspond to explosions of stars from relatively modest 
progenitor radii, and thus with intrinsic colours that are relatively 
red prior to peak. 

In the future, we will study explosions of WR stars with a 
more extended radius at death. For these, we expect bluer colours 
prior to maximum, and consequently different correction factors. 
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Table Al. Corrections factors for the 3p65Axl model sequence covering the photospheric phase. See Section[9]for discussion, as well as[Dessart & Hilliei 
<20051) for a description of the tabulated quantities and a presentation of the method used to calculate them. 


Age 

Tbv 

£bv 

Tbvi 

£bvi 

Tvi 

€vi 

T’phot 

-^phot 

^phot 

[day] 

[K] 


[K] 


[K] 


[K] 

[cm] 

[km s -1 ] 

2.360 

7760.0 

0.961 

11170.0 

0.505 

18040.0 

0.326 

7426.022 

3.863(14) 

18943.610 

2.600 

6460.0 

1.258 

9590.0 

0.561 

16330.0 

0.326 

6990.782 

4.064(14) 

18093.070 

2.860 

5600.0 

1.637 

8270.0 

0.664 

13940.0 

0.361 

6638.097 

4.216(14) 

17060.560 

3.150 

5040.0 

2.072 

7400.0 

0.781 

12310.0 

0.408 

6493.325 

4.328(14) 

15901.340 

3.470 

4680.0 

2.528 

6870.0 

0.899 

11360.0 

0.454 

6570.802 

4.373(14) 

14586.0 

3.820 

4420.0 

3.015 

6500.0 

1.017 

10730.0 

0.500 

6711.424 

4.388(14) 

13295.980 

4.200 

4250.0 

3.424 

6280.0 

1.104 

10400.0 

0.530 

6837.968 

4.441(14) 

12236.960 

4.620 

4130.0 

3.783 

6100.0 

1.191 

10120.0 

0.561 

6948.506 

4.496(14) 

11263.120 

5.080 

4060.0 

4.049 

5960.0 

1.273 

9760.0 

0.602 

7072.313 

4.535(14) 

10331.840 

5.590 

3960.0 

4.264 

5770.0 

1.340 

9310.0 

0.633 

7014.261 

4.625(14) 

9575.254 

6.150 

3860.0 

4.474 

5590.0 

1.396 

8910.0 

0.659 

6944.160 

4.765(14) 

8967.916 

6.760 

3810.0 

4.500 

5480.0 

1.417 

8620.0 

0.674 

6863.244 

4.981(14) 

8527.549 

7.440 

3790.0 

4.485 

5440.0 

1.411 

8510.0 

0.674 

6801.329 

5.268(14) 

8194.578 

8.180 

3820.0 

4.351 

5480.0 

1.381 

8550.0 

0.664 

6780.774 

5.615(14) 

7944.819 

9.0 

3920.0 

4.044 

5610.0 

1.319 

8700.0 

0.648 

6804.780 

6.029(14) 

7753.137 

9.900 

4050.0 

3.727 

5810.0 

1.242 

9060.0 

0.618 

6865.623 

6.505(14) 

7604.410 

10.890 

4230.0 

3.358 

6080.0 

1.155 

9550.0 

0.582 

6956.454 

7.035(14) 

7477.017 

11.980 

4450.0 

2.974 

6430.0 

1.053 

10170.0 

0.541 

7055.816 

7.629(14) 

7370.751 

13.180 

4720.0 

2.574 

6820.0 

0.956 

10970.0 

0.495 

7154.051 

8.282(14) 

7273.257 

14.500 

5030.0 

2.200 

7290.0 

0.853 

11710.0 

0.459 

7241.006 

8.974(14) 

7163.545 

15.950 

5340.0 

1.867 

7750.0 

0.756 

12610.0 

0.413 

7128.274 

9.891(15) 

7177.419 

17.550 

5640.0 

1.662 

8130.0 

0.710 

13200.0 

0.397 

7332.842 

1.044(15) 

6882.652 

19.310 

5870.0 

1.499 

8420.0 

0.664 

13520.0 

0.382 

7314.003 

1.118(15) 

6699.158 

21.240 

6010.0 

1.391 

8510.0 

0.643 

13270.0 

0.382 

7248.240 

1.189(15) 

6478.925 

23.360 

6010.0 

1.355 

8390.0 

0.643 

12680.0 

0.392 

7136.445 

1.256(15) 

6222.907 

25.700 

5750.0 

1.432 

7890.0 

0.684 

11490.0 

0.423 

6985.835 

1.323(15) 

5959.559 

28.270 

5450.0 

1.545 

7350.0 

0.740 

10350.0 

0.464 

6795.695 

1.376(15) 

5631.556 

31.100 

5120.0 

1.693 

6730.0 

0.828 

9070.0 

0.530 

6550.404 

1.400(15) 

5208.424 
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Table A2. Corrections factors for the 5pllAxl model sequence covering the photospheric phase. See Section [9] for discussion, as well as [Dessart & Hilliei 
<20051) for a description of the tabulated quantities and a presentation of the method used to calculate them. 


Age 

Tbv 

£bv 

Tbvi 

S.BVI 

T V i 

£vi 

Tphot 

^phot 

^phot 

[day] 

[K] 


[K] 


[K] 


[K] 

[cm] 

[km s -1 ] 

3.27 

3820.0 

2.660 

5590.0 

0.799 

9070.0 

0.369 

4895.1 

3.593(14) 

12717.7 

3.60 

3710.0 

3.040 

5340.0 

0.931 

8340.0 

0.442 

5322.9 

3.462(14) 

11129.6 

3.96 

3600.0 

3.227 

5150.0 

0.982 

7900.0 

0.470 

5245.4 

3.628(14) 

10603.8 

4.36 

3520.0 

3.383 

4990.0 

1.037 

7560.0 

0.497 

5267.4 

3.786(14) 

10050.8 

4.80 

3460.0 

3.428 

4900.0 

1.037 

7350.0 

0.502 

5085.7 

4.072(14) 

9817.7 

5.28 

3430.0 

3.543 

4810.0 

1.096 

7070.0 

0.543 

5226.7 

4.177(14) 

9155.7 

5.81 

3390.0 

3.680 

4760.0 

1.123 

7040.0 

0.547 

5218.4 

4.411(14) 

8786.7 

6.39 

3380.0 

3.758 

4780.0 

1.119 

7110.0 

0.543 

5221.0 

4.703(14) 

8517.7 

7.03 

3420.0 

3.625 

4840.0 

1.087 

7290.0 

0.520 

5122.0 

5.135(14) 

8453.9 

7.73 

3480.0 

3.584 

4960.0 

1.069 

7530.0 

0.511 

5262.9 

5.453(14) 

8164.8 

8.50 

3570.0 

3.314 

5100.0 

1.005 

7810.0 

0.479 

5111.8 

6.092(14) 

8295.6 

9.35 

3700.0 

3.140 

5290.0 

0.977 

8170.0 

0.470 

5340.3 

6.448(14) 

7981.9 

10.29 

3830.0 

2.912 

5500.0 

0.922 

8510.0 

0.451 

5390.7 

7.034(14) 

7911.9 

11.32 

3980.0 

2.591 

5720.0 

0.849 

8910.0 

0.419 

5281.0 

7.866(14) 

8042.9 

12.45 

4130.0 

2.441 

5910.0 

0.831 

9240.0 

0.415 

5503.4 

8.357(14) 

7769.5 

13.70 

4280.0 

2.244 

6120.0 

0.790 

9570.0 

0.401 

5565.2 

9.084(14) 

7674.5 

15.07 

4430.0 

2.079 

6350.0 

0.748 

9930.0 

0.387 

5622.4 

9.826(14) 

7546.9 

16.58 

4590.0 

1.915 

6550.0 

0.716 

10200.0 

0.378 

5685.1 

1.062(15) 

7414.2 

18.24 

4730.0 

1.796 

6740.0 

0.689 

10500.0 

0.369 

5756.9 

1.142(15) 

7247.7 

20.06 

4860.0 

1.700 

6930.0 

0.666 

10710.0 

0.364 

5827.7 

1.224(15) 

7062.5 

22.07 

5000.0 

1.599 

7090.0 

0.648 

10810.0 

0.364 

5893.5 

1.308(15) 

6859.3 

24.28 

5110.0 

1.535 

7190.0 

0.643 

10930.0 

0.364 

5971.6 

1.388(15) 

6614.2 

26.71 

5200.0 

1.480 

7230.0 

0.643 

10790.0 

0.373 

6018.1 

1.467(15) 

6358.7 

29.38 

5120.0 

1.521 

7000.0 

0.680 

10100.0 

0.405 

6010.5 

1.541(15) 

6072.5 

32.32 

5040.0 

1.567 

6790.0 

0.716 

9520.0 

0.438 

5996.6 

1.600(15) 

5729.8 

35.55 

4900.0 

1.654 

6530.0 

0.762 

8940.0 

0.474 

5962.6 

1.637(15) 

5330.4 

39.10 

4720.0 

1.750 

6170.0 

0.826 

8250.0 

0.520 

5853.2 

1.662(15) 

4918.9 

43.01 

4520.0 

1.841 

5800.0 

0.890 

7520.0 

0.575 

5695.2 

1.676(15) 

4510.5 

47.31 

4350.0 

1.855 

5440.0 

0.941 

6820.0 

0.630 

5539.1 

1.696(15) 

4148.1 

52.04 

4210.0 

1.796 

5120.0 

0.973 

6170.0 

0.684 

5386.5 

1.720(15) 

3826.0 

57.24 

4210.0 

1.590 

4950.0 

0.950 

5730.0 

0.716 

5282.2 

1.759(15) 

3556.5 

62.96 

4300.0 

1.343 

4890.0 

0.890 

5450.0 

0.721 

5233.8 

1.819(15) 

3344.1 

69.26 

4480.0 

1.091 

4910.0 

0.822 

5290.0 

0.707 

5212.4 

1.887(15) 

3153.6 
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Table A3. Corrections factors for the 6p5Axl model sequence covering the photospheric phase. See Section [9] for discussion, as well as[Dessart & Hillier 
J2005I) for a description of the tabulated quantities and a presentation of the method used to calculate them. 


Age 

Tbv 

£bv 

Tbvi 

iBVi 

Tvi 

£vi 

Tphot 

^phot 

^phot 

[day] 

[K] 


[K] 


[K] 


[K] 

[cm] 

[km s -1 ] 

3.89 

4140.0 

3.865 

5980.0 

1.288 

9510.0 

0.634 

7456.4 

2.991(14) 

8899.9 

4.28 

4030.0 

3.982 

5740.0 

1.347 

8900.0 

0.670 

7297.6 

3.126(14) 

8452.5 

4.71 

3880.0 

3.995 

5330.0 

1.442 

8170.0 

0.702 

6731.3 

3.317(14) 

8152.2 

5.18 

3790.0 

4.000 

5110.0 

1.496 

7790.0 

0.720 

6531.2 

3.518(14) 

7860.7 

5.70 

3740.0 

4.000 

5000.0 

1.514 

7600.0 

0.724 

6403.5 

3.762(14) 

7639.7 

6.27 

3710.0 

3.995 

4960.0 

1.505 

7490.0 

0.724 

6318.4 

4.036(14) 

7450.4 

6.90 

3690.0 

4.000 

4950.0 

1.491 

7420.0 

0.724 

6265.5 

4.337(14) 

7274.3 

7.59 

3690.0 

4.000 

4980.0 

1.464 

7450.0 

0.715 

6242.4 

4.660(14) 

7105.5 

8.35 

3700.0 

4.000 

5030.0 

1.437 

7490.0 

0.711 

6244.4 

5.009(14) 

6943.2 

9.19 

3720.0 

3.995 

5100.0 

1.405 

7580.0 

0.702 

6268.5 

5.390(14) 

6788.7 

10.11 

3740.0 

4.000 

5200.0 

1.365 

7690.0 

0.693 

6306.2 

5.808(14) 

6649.1 

11.12 

3800.0 

3.896 

5310.0 

1.324 

7850.0 

0.679 

6359.4 

6.245(14) 

6500.4 

12.23 

3890.0 

3.689 

5420.0 

1.288 

8010.0 

0.665 

6410.8 

6.726(14) 

6365.4 

13.45 

4000.0 

3.450 

5560.0 

1.239 

8150.0 

0.656 

6471.9 

7.213(14) 

6207.1 

14.80 

4100.0 

3.251 

5690.0 

1.193 

8320.0 

0.643 

6524.5 

7.733(14) 

6047.1 

16.28 

4200.0 

3.057 

5800.0 

1.157 

8420.0 

0.634 

6550.7 

8.278(14) 

5885.1 

17.91 

4300.0 

2.867 

5910.0 

1.121 

8510.0 

0.625 

6576.9 

8.860(14) 

5725.8 

19.70 

4390.0 

2.700 

5980.0 

1.094 

8540.0 

0.620 

6585.9 

9.458(14) 

5556.5 

21.67 

4460.0 

2.574 

6040.0 

1.072 

8520.0 

0.620 

6585.2 

1.008(15) 

5381.5 

23.84 

4520.0 

2.470 

6090.0 

1.054 

8470.0 

0.625 

6574.4 

1.069(15) 

5192.2 

26.22 

4580.0 

2.371 

6100.0 

1.049 

8370.0 

0.634 

6555.7 

1.130(15) 

4989.7 

28.84 

4610.0 

2.312 

6090.0 

1.049 

8230.0 

0.647 

6522.3 

1.193(15) 

4789.5 

31.72 

4620.0 

2.276 

6070.0 

1.049 

8120.0 

0.656 

6474.0 

1.256(15) 

4583.8 

34.89 

4630.0 

2.213 

6010.0 

1.054 

7900.0 

0.675 

6394.5 

1.324(15) 

4393.7 

38.38 

4490.0 

2.290 

5750.0 

1.108 

7420.0 

0.720 

6240.1 

1.397(15) 

4214.1 

42.22 

4400.0 

2.299 

5590.0 

1.126 

7150.0 

0.733 

6090.4 

1.472(15) 

4034.7 

46.44 

4240.0 

2.358 

5330.0 

1.162 

6660.0 

0.778 

5863.9 

1.550(15) 

3862.5 

51.08 

4100.0 

2.367 

5110.0 

1.175 

6320.0 

0.792 

5665.2 

1.622(15) 

3676.1 

56.19 

4020.0 

2.222 

4920.0 

1.153 

5950.0 

0.801 

5453.5 

1.699(15) 

3499.6 

61.81 

4020.0 

1.924 

4780.0 

1.090 

5580.0 

0.805 

5268.7 

1.778(15) 

3330.0 

67.99 

4090.0 

1.586 

4690.0 

1.013 

5260.0 

0.805 

5140.7 

1.855(15) 

3158.2 

74.79 

4230.0 

1.270 

4660.0 

0.927 

5040.0 

0.792 

5063.5 

1.933(15) 

2992.2 
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